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1) Idée sur la méthode 


Ce projet-là présente une méthodologie d’optimisation de l’enveloppe d’un 
bâtiment du point de vue de la performance énergétique. Des critères permettant 
d'évaluer la performance énergétique sont définis, comme la consommation 
énergétique annuelle du bâtiment et le confort thermique des habitants. 
L’optimisation est réalisée grâce à un couplage entre un réseau de neurones 
artificiels et un algorithme génétique. 


2) Notions de performance énergétique 


a) Introduction 


Cette méthode, présenter par le Grenelle de l’Environnement, propose la mise 
en œuvre à grande échelle d’un programme pour réduire les consommations 
énergétiques des bâtiments. Donc, on doit jouer sur l’enveloppe des bâtiments 
pour réaliser des constructions performantes. Vince et Monchoux ont réalisé 
l’optimisation d’une paroi multicouche seule d’un point de vue de sa résistance 
thermique et de sa capacité thermique. Dans ce travail-ci, nous souhaitons changer 
d’échelle et mener l’optimisation de l’enveloppe d’un bâtiment. Il convient 
d’intégrer l’occupant et de savoir si la nature des parois extérieures influe sur son 
confort. 


Alors on va chercher d’optimiser la paroi opaque verticale extérieure d’un 
bâtiment permettant de le rendre performant d’un point de vue énergétique. Nous 
allons chercher les propriétés thermo physiques équivalentes du mur extérieur 
(conductivité thermique K, chaleur spécifique volumique pc) qui maximisent la 
performance énergétique. La notion de performance énergétique doit combiner à 
la fois l’aspect énergétique et le confort de l’occupant. Ainsi, nous devons 
déterminer des indicateurs appropriés qui traduisent cette notion. 


b) Revue bibliographique 


Une étude d’optimisation a besoin de fonctions objectifs (FO) qu’il faut 
définir. Ces fonctions objectives doivent être minimisées sur l’espace de 
recherche. Dans le cadre de ce travail, elles doivent caractériser précisément la 
performance énergétique du bâtiment. Il en existe un certain nombre qui sont 
proposés dans la littérature. 


Une étude propose la caractérisation du confort thermique par le degré 
d’inconfort en été et en hiver Lise et lhiver. Ces 2 grandeurs sont définies par : 


kte = [one (Tin — Th)dt Tin > Ty 
lhiver = fagar = Tin)dt TL > Tin 


Où : 

Tin: temperature interieur 

Ty: température de confort d’été 
T,: température de confort l’hiver 


Iste et Liver Comptabilisent l’accumulation des degrés excédant la 
température de confort. 
Compte tenu des définitions, le confort est d’autant plus assuré que Iere et 
Thiver tendent vers 0. 


Wright et Nancy réalisent une optimisation qui porte simultanément sur le coût 
énergétique d'un système CVC (chauffage, ventilation et climatisation), et le 
minimum du PPD (Pourcentage prévu d'insatisfaits). Leur optimisation est 
effectuée par un algorithme génétique. Magnier et Haghighat cherchent également 
à optimiser à la fois la consommation énergétique et le confort des occupants. La 
fonction-objectif qu’ils choisissent pour le confort thermique est basée sur le 
PMV (Vote moyen prévu). 


Dans la méthodologie d’optimisation que nous proposons, nous allons 
également utiliser les algorithmes génétiques. De plus, il nous semble 
indispensable que dans la conception de performance énergétique du bâtiment 
soient associées sa consommation énergétique et le confort d’été. Ce dernier sera 


caractérisé par le degré d’inconfort Isté 
Justification du choix 


Etant futurs ingénieurs en génie civil, on est sensé d’améliorer les conditions de 
vie des citoyens en tenant compte bien sûr de l’énergie, qui devient de plus en plus 
rare, et de l’environnement, qui ne peut plus supporter les émissions. Alors, on 
doit trouver une solution qui va en même temps satisfaire les besoins des 
consommateurs sans nuire à l’environnement pour garantir notre vie et celles des 
générations futures. Donc, c’est là ou commence l’idée de cette présentation afin 


d’atteindre le but de garantir un confort pour les consommateurs sans nuire à 
l’environnement. 


3) Méthodologie d'optimisation de la performance énergétique 


La restriction majeure des algorithmes génétiques, et particulièrement des 
algorithmes génétiques multi objectifs, est leur nécessité de centaines voire de 
milliers d'évaluations pour atteindre une solution optimale. Dans le domaine du 
bâtiment, où les évaluations sont généralement réalisées par des simulations 
thermiques dynamiques relativement chronophages (TRNSYS, EnergyPlus, 
Comfie...), le temps nécessaire à l’optimisation peut s’avérer élevé et devenir 
rapidement trop cher 


Une méthode efficace pour réduire considérablement ce temps, est de réaliser 
une approximation par surface de réponse. La surface de réponse est une fonction 
liant les réponses du système aux variables d’entrée. Ainsi, l’approximation par 
surface de réponse consiste à substituer à la vraie surface, une surface définie 
analytiquement, plus facile à optimiser. La détermination de cette surface de 


réponse s'effectue au moyen d’un réseau de neurones artificiels multicouche. 


a) Réseau de neurones artificiels 


Un réseau de neurones artificiels (RNA) est composé d’éléments simples, des 
neurones, fonctionnant en parallèle. Son principe a été fortement inspiré par le 
système nerveux biologique. Les informations d’entrée circulent à travers 
plusieurs couches de neurones afin de délivrer au final un signal de sortie. Plus 
précisément, un RNA est composé d’une couche de neurones d’entrée (variables 
du système), d’une couche de sortie (réponse approximée) et d’une ou plusieurs 
couches cachées reliant les couches d’entrée et de sortie. Le fonctionnement du 


RNA est fortement influencé par la connexion des neurones entre eux 








Couche d’entrée couche cachée couche de sortie 


Le RNA doit être entraîné à effectuer une tâche spécifique en ajustant les poids 
entre les neurones. Cet apprentissage des RNA est réalisé de sorte que, pour une 
entrée particulière présentée au RNA, corresponde une cible spécifique. Les poids 
sont ajustés par comparaison entre la réponse du RNA et la cible, jusqu’à ce que la 
sortie corresponde au mieux à la cible. Ce type d’apprentissage est dit supervisé et 
nécessite un nombre relativement élevé de paires « Entrée/sortie ». 


b) Algorithmes génétiques NSGA II 


Les algorithmes génétiques sont des algorithmes d'optimisation s'appuyant 
sur des techniques dérivées de la génétique et des mécanismes d'évolution de la 
nature : au sein d’une population, seuls les individus les mieux adaptés ont une 
descendance à qui ils transmettent leurs gènes. L’espèce s’améliore ainsi par 
croisement au fil des générations, et aussi grâce au phénomène de mutation des 
gènes. Depuis le milieu des années 1990, les algorithmes génétiques ont fait leur 
preuve dans le domaine de la thermique. 


Nous avons choisi dans ce travail d'utiliser l’algorithme génétique multi 
objectif NSGA II (Non-dominated Sorting Genetic Algorithm). L’approche 
NSGA II est basée sur le classement non dominé des individus (notion de 
dominance de Pareto) et comporte des stratégies élitistes pour conserver les 
meilleures solutions trouvées au cours de l’exécution de l’algorithme. 


c) Structure de méthodologie 


La méthodologie d'optimisation proposée est schématisée sur la Figure ci- 
dessous. D’abord, un modèle de bâtiment est réalisé sous TRNSYS (un logiciel de 
simulation thermique dynamique des systèmes). Ensuite, à partir de ce modèle, est 
générée une base de données qui va servir à l’apprentissage du RNA. Cette base 
de données est constituée des valeurs des variables d’optimisation, à savoir les 
propriétés thermo physiques de l’enveloppe (K et pc) et des deux fonctions- 
objectifs correspondantes calculées par TRNSYS, c’est-à-dire la consommation 
énergétique annuelle notée Qror et le degré d’inconfort Fete. Elle est obtenue par 
le logiciel GenOpt qui permet la création des données suffisamment réparties sur 
l’espace de recherche. Une fois l’entraînement et la validation assurés, le RNA est 
prêt à effectuer des évaluations rapides des fonctions-objectifs avec une bonne 
précision, sans dégrader la complexité du modèle. Enfin, NSGA II est couplé 
avec le RNA afin de réaliser l’optimisation. 


En effet, NSGA II fournit des valeurs d’entrée au RNA et le RNA effectue 
les évaluations dont a besoin le NSGA II. 






Simulations (TRNSYS 
+ GenOpt) 






RNA validé 


Fonctions-objectifs 


NSGA II 













Apprentissage et 
validation du RNA 


4) Application 













Solutions optimales 






Nous allons appliquer la méthodologie développée à un cas de maison 
individuelle, dite Maison Mozart. Nous ne donnons ici que les principaux 
éléments utiles à cette étude. Il s’agit d’une maison de plain-pied d’une surface 
hors œuvre nette de 112 m°. La hauteur du plafond est de 2.3 m. Sur le plan ci- 
dessous, les murs en noir correspondent aux murs extérieurs. 




















Les zones thermiques ont été modélisées par secteur d’activité similaire et 
orientation des pièces. La surventilation nocturne en été est intégrée au modèle : 
elle est déclenchée la nuit lorsque la température d’air intérieure est supérieure à 
celle de l’air extérieur. 


La paroi extérieure est composée d’une couche massive de 20 cm et d’une 
couche fixe d’isolant de 8 cm positionnée à l’extérieur. 


Conformément à notre objectif, ce sont les propriétés thermo physiques (k, 
pc) de la couche massive que l’on cherche à optimiser du point de vue de la 
consommation énergétique annuelle et du confort d’été. La consommation 
énergétique annuelle Qror est le besoin de chauffage pour maintenir la 
température de consigne de chauffage (19°C) en hiver. Nous rappelons que 
l’inconfort d’été est caractérisé par Ięté. La température de confort d’été Ty est 


fixée à 28°C. Ces 2 grandeurs sont calculées par TRNSYS. 
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Nous allons faire varier k et pc entre les bornes données par le Tableau 1. Le 
choix de ces bornes est basé sur l’observation des matériaux usuels de 
construction. La limite basse de k se justifie par la présence de l’isolant dans la 
paroi extérieure. Ainsi, la couche massive n’a pas à intégrer la totalité de la 
fonction isolante de la paroi. 
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min 0.05 4x10* 





max 1.75 2x106 

















Tableau 1 : Gammes de variation de k et pc 


De plus, nous faisons également varier le climat afin de noter son influence 
sur l’optimisation. Nous optons pour 2 climats : un climat chaud, Nice et un climat 
froid, Nancy. 


a) Entraînement et validation du RNA 


Le RNA comporte 2 neurones dans la couche d’entrée (k, pc), 20 neurones 
dans la couche cachée et 2 neurones dans la couche de sortie (Qror, Lete). La phase 
d’entraînement du RNA nécessite 130 simulations réparties de façon aléatoire par 
GenOpt. Ces simulations fournissent Qror et ls de maisons Mozart constituées 
de paroi extérieures différentes. Une fois l’entraînement réalisé, nous vérifions que 


le RNA concorde bien avec la modélisation de TRNSYS. En se basant sur un 
échantillon de 10 combinaisons de paramètres d’entrée, nous obtenons un écart 
relatif moyen entre la valeur de Qror du RNA et de TRNSYS de 1.2 %. Pour Lite, 
l’écart relatif moyen est de 4.6 %. Nous en concluons que notre RNA est validé et 
qu’il est donc prêt à être couplé à l’algorithme génétique NSGA IT. 


b) Optimisations et résultats 


Une étude de sensibilité numérique a été réalisée au préalable pour déterminer 
les paramètres de l’algorithme génétique NSGA II. Les individus sont les 
différentes parois extérieures caractérisées par leurs propriétés thermo physiques 
k et pc, nous en avons choisi 50. Il est nécessaire de suivre l’évolution de 
l’algorithme sur 300 générations pour atteindre la solution optimale convergée. 


Les figures ci-dessous présentent les résultats de l’optimisation pour le climat de 
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Nous constatons que les fronts de Pareto obtenus sont bien étalés. De plus, 
pour les 2 climats, ils ont un comportement similaire, même si les valeurs sont 
très différentes. Nous constatons également 2 zones distinctes : une zone 
délimitée par les points A et B pour Nancy (respectivement, A’ et B’ pour Nice) 
et une zone entre B et C (respectivement, B’ et C’). 


La 1° zone (de A à B et de A’ à B’) correspond à une brusque croissance de 
Le pour les plus faibles valeurs de Qror. En effet, ce sont les parois extérieures 
les plus isolées (valeurs minimales de k) qui sont les solutions optimales et de fait, 
le confort d’été est nettement dégradé, car la chaleur ne peut suffisamment être 


évacuée à l’extérieur. De plus, nous observons que, dans cette zone, la variation 
de pc a une influence négligeable sur Qror. 



































0.05 893.0 1414 
Nice B’ 0.05 1987.2 1445 346 
C 1.75 2000.0 1885 95 





Tableau 2 : Valeurs associées aux points À, B, C, A’, B’ et C’ 


La 2°% zone (de B à C et de B’ à C’) est beaucoup plus large et est quasiment 
linéaire. Nous constatons que pc est globalement constant pour toutes les 
solutions de cette zone et est fixé à sa valeur maximale, 2000 kJ.m?.K!. 


Beaucoup de solutions optimales obtenues ne correspondent pas à des 
matériaux existants : les valeurs de k sont plus élevées pour les matériaux réels. Si 
nous désirons des solutions plus réalistes, il faudrait réaliser l’optimisation avec 
une base de données de matériaux usuels de construction. 


5) Conclusion 


Ce travail propose une méthodologie d’optimisation de l’enveloppe d’un 
bâtiment du point de vue de sa performance énergétique. Cette méthodologie 
nécessite en amont une recherche précise de la caractérisation des critères de 
performance énergétique. Afin de réduire le temps de calcul sans ôter la 
complexité du problème, un réseau de neurones artificiels a été élaboré : il doit 
être au préalable entraîné et validé. Son rôle consiste à fournir des évaluations 
précises et rapides des fonctions-objectifs à un algorithme génétique de type 
NSGA II. 

Cette méthodologie est testée pour une maison individuelle où il s’agit de 
trouver les propriétés thermo physiques optimales de la paroi extérieure qui 
minimisent à la fois Qror et Tete pour Nancy et Nice. 


6) Application au Maroc 


Le Maroc, comme tous les autres pays du monde, il a des problèmes énergétiques 
alors cette méthode sera bénéfique pour lui or il existe un grand problème, tout d’abord 
celle de l’argent parce qu'une telle étude demande une somme importante d’argent, 
chose qui n’est pas offerte à tout le monde, ensuite 1l y a la contrainte du temps qui joue 


un rôle primordial pour le cout et enfin les matériaux utilisés qui sont rares. 


Pour conclure, on peut dire que cette méthode est applicable au Maroc sauf qu’il 


faut essaye de l’optimiser plus afin de la rendre accessible à tout le monde en jouant sur 
la dure de l’étude et de trouver des matériaux qui ont les mêmes exigences demander 
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